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. I. INTRODUCTION 
(N ; 

' Gemini surfactants arc composed of two monomeric surfactant molecules linked by a spacer chain. They constitute 

] a new class of amphiphilic molecules having its own distinct behavior. Since their first systematic studies over a decade 

'■^ ■ ago, gemini surfactants have been the subject of intensive research (see Ref. JJ and references therein). Research has 

' been motivated by the advantages of gemini surfactants over regular ones with respect to various applications, e.g., 

. their increased surface activity, lower critical micelle concentration (cmc), and useful viscoelastic properties such as 

^ effective thickening. 

Besides their importance for applications, the behavior of gemini surfactants is qualitatively different in several 

Q ■ respects from that of regular surfactants, posing challenges to current theories of surfactant self-assembly. The main 

I puzzles raised by gemini surfactants can be summarized as follows Q. 

, • Surface behavior. The area per molecule in a saturated monolayer at the water-air interface, made of gemini 

■ surfactants with polymethylene spacers [m-s-m surfactants, where s is the spacer length and m the tail length 

I ' in hydrocarbon groups), has a non-monotonous dependence on s For example, for tail length of to = 12 

''O ■ the molecular area at the water-air interface is found to increase with s for short spacers, reach a maximum 

Ch ' at about s ~ 10-12, and then decrease for longer spacers. This decrease in the specific area for the m-s-m 

O I surfactants is somewhat unexpected given the fact that the molecule becomes bigger as s increases. One would 

1 I ' naively expect a monotonous increase in the molecular area as indeed is observed for another class of gemini 
I surfactants having a poly (ethylene oxide) spacer (to-EOjj-to surfactants) 0,0]. 

I ! 

^ , • Micellization point. The critical micelle concentration (cmc) of gemini surfactants is typically one to two orders 

' of magnitude lower than that of the corresponding monomeric surfactants having the same head and tail groups 
For regular (monomeric) surfactants the cmc decreases monotonously with the number of hydrocarbon 

• groups because of increased molecular hydrophobicity. In the case of m-s-m gemini surfactants, by contrast, the 

I dependence of the cmc on the spacer length s is non- monotonous with a maximum at about s ~ 4-6 H . |8[ . 

. Similarly, the Krafft temperature exhibits a minimum |^ and the micellization enthalpy a maximum |lO| at 

' about the same s value. 

> ' 

^ ; 

S ' of decreasing curvature, e.g., from spherical micelles to cylindrical micelles to bilayer vesicles |llllT^. However 

. when the polymethylene spacer length in m-s-m gemini surfactants is increased, a different sequence of shapes 

' is observed, for instance, from cylindrical micelles to spherical micelles to vesicles for the 12-S-12 surfactants 

O , |l3lll4j . Moreover, gemini surfactants with short spacers exhibit uncommon aggregate morphologies in the form 

O ' of branched cylindrical micelles and ring micelles \la^ . 

. ^ ' • Phase behavior. The spacer length in m-s-m gemini surfactants has an unusual effect also on the phase behavior 

. of binary surfactant -water mixtures. For geminis with tail length to = 12, for instance, the phase-diagram region 

^ ' corresponding to hexagonal and lamellar phases is found to shrink with increasing s, disappear for s = 10-12, 

5^ , and re-appear for s = 16 p^ . In ternary systems of water-oil-TO-s-TO surfactant, the size of the microemulsion 
|smgle-phase) region in the phase diagram has a non-monotonous dependence on s with a maximum at s ~ 10 

• Dynamics. Dilute micellar solutions of gemini surfactants with short spacers have unusual rheological properties, 
such as pronounced increase in viscosity upon increase of surfactant volume fraction p^ll9| | and shear-thickening 

MM. 



• Aggregate shape. As certain parameters, such as the relative size of the head and tail groups or the salt 
concentration, are progressively changed, regular surfactants change their aggregate morpholo gy i n the direction 
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In view of the amount of experimental work and its unusual findings, the number of theoretical studies devoted to 
gemini surfactants has been surprisingly small. In this chapter we have, therefore, two aims. The first is to review 
the current state of theoretical models of gemini surfactants. The second, perhaps more important aim, is to indicate 
the considerable gaps in our knowledge and the key open questions awaiting theoretical work. In Sec.^]wc set the 
stage by reviewing several theoretical models of surfactant self-assembly. This will facilitate the discussion in Sec. IIIII 
of the gemini surfactant models, which can be viewed as extensions to previous models of regular surfactants. Finally, 
in Sec. II VI we conclude and summarize the open questions. 



II. MODELS OF SURFACTANT SELF-ASSEMBLY 



In this section we review several theoretical models pertaining to the self-assembly of soluble surfactants. This is 
not meant to be an exhaustive review of self-assembly theory but merely to provide the appropriate background for 
the models of gemini surfactants discussed in Sec. IIIII 



A. Surface behavior 



Let us start by considering an aqueous surfactant solution below the cmc. The soluble surfactant molecules self- 
assemble into a condensed layer at the water-air interface, referred to as a Gibbs monolayer (to be distinguished 
from Langmuir monolayers that form when insoluble surfactants are spread on the water-air interface) p3 |. Since 
the surfactant is water-soluble, this layer exchanges molecules with the bulk solution and a nonuniform concentration 
profile forms. Typical surfactants have strong surface activity, i.e., the energy gained by a molecule when it migrates 
to the surface is much larger than the thermal energy fceT. As a result, the concentration profile drops sharply to its 
bulk value within a molecular distance from the surface (hence the term monolayer). 

The number of molecules participating in a Gibbs monolayer per unit area, the surface excess F, is obtained from 
integrating the excess concentration (with respect to the bulk one) over the entire solution. Such a monolayer can be 
regarded as a separate sub-system at thermodynamic equilibrium and in contact with a large reservoir of molecules 
having temperature T and chemical potential /i. From the excess free energy per unit area of this system, ^(T,ii), 
which is by definition the surface tension of the solution, we get the number of molecules per unit area: 




This is referred to as the Gibbs equation (2^. For dilute solutions fj. oc fceTlnC, where C is the bulk surfactant con- 
centration. (The constant of proportionality is one for nonionic surfactants and ionic ones at high salt concentration; 
it has a higher value for salt-free ionic surfactant solutions, where strong correlations between the different ions lead 
to non- ideal activity coefficients .22].) Hence, 

Because of the high surface activity of surfactant molecules, leading to a sharp concentration profile at the water-air 
interface, F"^ is commonly interpreted as the average surface area per molecule, a. The second consequence of the 
high surface activity is that, already for C much smaller than the cmc, the monolayer becomes saturated, i.e., F stops 
increasing with C. Experimentally, the curve describing the change in 7 as a function of In C becomes a straight line 
with a negative slope proportional to — F. 

We now wish to find a simple estimate for the energy of lateral interaction between molecules in such a saturated 
monolayer (repeating a well-known result of Ref. [llj). Saturation implies that the molecules are packed in an 
energetically optimal density, such that there is no gain in adding or removing molecules. This optimum arises from 
a competition between attractive and repulsive interactions. The attractive interaction tries to decrease the area 
per molecule, and we can phenomenologically write its energy per molecule as proportional to a, 71a, where the 
proportionality constant 71 has units of energy per unit area. Since the attraction comes mainly from the desire of 
the hydrocarbon tails to minimize their contact with water, 71 should be roughly equal to the hydrocarbon-water 
intcrfacial tension (71 ~ 50mN/m). The repulsive interaction, on the other hand, tries to increase a and, at the 
same phenomenological level, we can write its energy per molecule as inversely proportional to a, a/a, where a is 
a positive constant. Minimizing the sum of these two contributions we get for the interaction energy per molecule 
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u = "/i{a — ao)^/a + const, where oq = is the optimal molecular area. Expanding around a = oq to second 

order and omitting the constant term, we obtain 

u{a) ~ —(a - ao)^. (3) 

ao 

This is merely a harmonic approximation for the interaction energy associated with small deviations from optimal 
packing. 

We can slightly modify this result to obtain a similar harmonic estimate for the energy U2 of effective interaction 
between two neighboring molecules residing in the saturated monolayer. (This will be useful later on when we add the 
spacer to form gemini surfactants.) We need to divide u by half the number of inplane neighbors, q/2, and express a 
in terms of the average intermolecular distance a — rjr^ , where 77 is a prefactor of order unity (e.g., for hexagonal 
packing q — & and rj — V3/2 ~ 0.9). Assuming again that r is close to its optimal value ro, we get 

U2(r)^ifco(r-ro)^ A^o = (4) 
2 q 

This expression replaces the actual surfactant-surfactant interaction in the monolayer, including effects of other nearby 
surfactants, with an "effective spring" of equilibrium length rg and spring constant fco. For hexagonal packing we get 
the reasonable value ~ 120mN/m ~ 0.3 /cbT/A^. Note that, because of the saturation condition, the expression 
for fco is insensitive to molecular details. In turn, those details will affect the properties of the saturation state itself, 
e.g., the value of oq or tq. 



B. Micelles 



As the solution concentration is increased beyond the cmc the surfactant molecules start to form aggregates. Unlike 
simple solute molecules (e.g., alkanes), which undergo macroscopic phase separation upon increasing concentration 
or changing temperature, surfactants form micelles at the mesoscopic scale. The challenges posed to theories of 
surfactant self-assembly are to predict the micellization point as a function of concentration (i.e., the cmc, hereafter 
referred to by the corresponding volume fraction (pcmc) and temperature (Tin), as well as the micelle shape and size. 
The main complications come from the fact that micellization is not a macroscopic phase transition — the aggregate 
sizes are finite and polydisperse — and thus the well-developed theoretical framework of phase transitions does not 
strictly apply. 

From a thermodynamic point of view, the difference between surfactant micellization and phase separation lies in 
the following observation [1^]. For alkanes solubilized in water, for example, the (Gibbs) free energy per molecule in 
an aggregate of size N , g^, is a monotonously decreasing function of — for A^ — > cx) gjv tends to the free energy 
per molecule in the bulk alkane phase, goo, while for smaller N g^ > goo due to unfavorable surface terms of the 
finite cluster. As a result, there is a critical concentration (or critical temperature) at which the favorable size changes 
discontinuously from monomers solubilized in water (A^ = 1) to a macroscopic phase of bulk alkane (A^ — + 00). The 
first-order phase transition point is reached when the chemical potential of monomers exceeds minj^Tv} = <?oo- In a 
dilute solution this implies that ipc = e^°°^'^''^'^\ fceTc = goo/^^f- (We have set the free energy of the A^ = 1 state as 
the reference, 171 = 0.) In the case of surfactants, by contrast, gN has a minimum at a finite aggregate size A^*. As a 
result, when the chemical potential exceeds gN', a large population of aggregates appears, whose sizes are distributed 
around A^*. Hence, the micellization point can be estimated as 

^cmc = e»"* /('^^^^ , fcBT„, = gN' I In ^. (5) 

The remaining task is to obtain a theoretical expression for g\j , the free energy per molecule in aggregates of size A^ 
and shape S. Since we expect this function to have a minimum at a finite yet large A^ (say, A^* ~ 10^), the importance 

(S) 

of many-body interactions is inevitable, and obtaining g\^ from rigorous statistical mechanics is a formidable task. 
Consequently, analytical models have relied on phenomenological approaches, trying to account for various competing 
contributions to the free energy while assuming a certain geometry for the aggregate (111 |12l i23l] . From the minimum 

(s) 

oig]^ with respect to A^ and various possible shapes S one can obtain the aggregate shape, aggregation number, and 
micellization point, using Eq. (O. 

In the simplest picture [TH Il2j| rough estimates for the minimum of g^^ can be obtained by imposing geometrical 
constraints that arise from the incompressibility of the micellar hydrocarbon core (see Fig.^. These constraints lead 
to a finite aggregation number A^*, as required. The hydrocarbon tail chains cannot extend beyond a certain length 
I, and each tail must occupy a certain volume v. (Both I and v are known to have a simple linear dependence on the 
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number of hydrocarbon groups in the tail chain [2j-) addition, as in the case of the Gibbs monolayer of Sec. Ill Al 
surfactant molecules in the favorable aggregate state attain a certain optimal area per molecule Oq. Then, due to 
translational entropy, the favorable aggregates would be the smallest ones that satisfy all constraints. The constraints 
define a dimensionless packing parameter (Fig. , 

P-^- (6) 

If P < 1/3 the constraints can be satisfied by spherical micelles, which will be the smallest and hence the most 
favorable ones; when l/3<P<l/2 the micelles will have to be of elongated or cylindrical shapes; for 1/2 < P < 1 
planar shapes will form; and for P > 1 the morphology must be inverted. 




FIG. 1: Packing constraints on a surfactant in an aggregate. Each head group occupies an optimal area ao on the aggregate 
surface; the tail chain occupies a volume v and cannot stretch beyond length I. These constraints define the packing parameter, 
P = v/{aol), which suggests the possible aggregate shape. 



Once the shape is determined we can find the maximum allowed aggregation number. For example, for spherical 
micelles N < AirP/ao — {4:Tr/3)P/v, and we get 

TV < 367r^. (7) 

We see how competing interactions between the molecules (giving rise to ao) together with the incompressibility of 

the micellar core lead to finite micelles. Since the tail chains usually should not stretch to their full extent, the actual 

aggregation number will be smaller than this upper bound. 

Yet, these geometrical arguments cannot provide us with theoretical predictions as to the optimal molecular area 

(s) 

ao itself or the aggregation free energy gjsfi , as well as their dependence on parameters such as temperature or salt 
concentration. In order to get such information and subsequently predict the micellization point, micelle shape and 
size, one needs a more detailed theory. 



1. Phenomenological models 

There have been attempts to analytically account for the various competing contributions to the free energy per 
molecule gj^^ (e.g., Ref. [l^]). The advantage of this approach is that, once we have an expression for the free 
energy, we can easily change parameters and gain insight into the role of various contributions. On the other hand, 
such models essentially attempt to push the limits of the phenomenological approach toward a detailed molecular 
description. They usually entail uncontrolled approximations and parameters whose accurate values are often hard 
to obtain. As an example, which will serve us in Sec. lIIII we give here an analysis along the lines of (yet not identical 
to) Ref. [23. 

Five major contributions to the aggregation free energy (per surfactant molecule on the aggregate) can be considered 



57V (") = 5hc + 5int(a) + .gcs(a, R) + gst(a) + 9c\{R), 



(8) 
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where a is the area per molecule on the aggregate surface and R the aggregate size (radius or width). Note that R 
is not an independent variable but is related to N and a via the aggregate geometry S, e.g., for spherical micelles, 
Na = AttR^. 

(i) The driving force for aggregation is the hydrophobic effect, i.e., the free energy per surfactant molecule ghc 
gained by shielding the hydrocarbon groups from water This contribution to 5^'' is negative and, to a good 
approximation, independent of N and the aggregate geometry S. Namely, its contribution to the entire aggregate 
free energy is linear in N and tends to increase the aggregate size. The hydrophobic term ghc depends linearly on 
the number of hydrocarbon groups in the surfactant, with a reduction of roughly k^T per hydrocarbon group |l2j |. 
That is why, for regular surfactants, the cmc decreases exponentially with the number of hydrocarbon groups in the 
molecule and is reduced by a factor of roughly 2-3 per each additional hydrocarbon group. 

(ii) The hydrophobic gain is corrected by an interfacial contribution g^^i due to the unfavorable contact between 
the hydrocarbon core and water. 



5int (a) = 7i (a - 



0, 



(9) 



where 71 is the interfacial tension of the core-water interface (roughly equal to the hydrocarbon-water interfacial 
tension), and flmin is the minimum area per molecule, i.e., the interfacial area occupied by a head group. This 
contribution evidently acts to reduce the area per molecule. 

(iii) If the surfactant head groups are charged, there is electrostatic repulsion between them, acting to increase a. 
Within the Poisson-Boltzmann theory this electrostatic contribution is given by : 



<?es(a, R) - 2fcBT { ln[/3 + (1 + /32)i/2] _ [(1 + ^2^1/2 _ 



2cAr 



In 



(10) 



where (3 = 47rZBAD/a is a dimensionless charging parameter depending on two other lengths, the Debye screening 
length Ad and the Bjerrum length Zb- The Debye screening length in the solution is Ad — (87rZBCsait)~^^^, where Csait 
is the added salt concentration, taken here to be monovalent, and Zb = / i^k-QT) is about 7 A for aqueous solution 
with dielectric constant e = 80 at room temperature. (For simplicity, a monovalent head group has been assumed.) 
Finally, c is the mean curvature of the aggregate (e.g., 1/R for spherical micelles). 

(iv) There is also steric repulsion between head groups. From the (non-ideal) entropy of mixing per molecule we 
get for this contribution 



gst{a) = fceT [ln(a,nin/a) + (a/amin - 1) ln(l - amin/a)] 



(11) 



(v) The last contribution to the free energy is associated with the tail packing in the hydrophobic core, i.e., deviations 
of the hydrocarbon tail chains from their relaxed length Iq, 



9el{R)^\k'{R~lof. 



(12) 



The "elastic constant" k' depends on the chain statistics, as well as the packing parameter (i.e., aggregate shape) [23. 

The equilibrium aggregation number N and specific area a (and hence also aggregate size R) for a given shape S 
and surfactant chemical potential fi are then determined by the equations 



da 



= 0. 



(13) 



Comparing the minimum value of g]^ for various shapes S, one also obtains the equilibrium aggregate morphology. As 
long as fj, < kBThiipcmc-, these equations will have no solution, and the monomeric state (iV = 1) of single surfactant 
molecules solubilized in water is the stable one. As the chemical potential increases, we reach the micellization 
condition given by Eq. jSJ, where the average micelle size at the cmc, N* , can be calculated now from the expression 

of at its minimum. 



2. Computer simulations 

Another route to overcome the complexity of treating surfactant micellization is to use computer simulations. This 
approach can be divided into two categories: statistical-mechanical models using Monte Carlo (MC) simulations, and 
Molecular Dynamics (MD) simulations. 
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Following Widom's statistical- mechanical model of microemulsions |2g , a host of lattice models were presented for 
treating surfactant self-assembly (see, e.g., Refs. [S^, IH, |33, |^ ) . These molecular "toy models" represent the 
water molecules and various groups in the surfactant as Ising spins on a discrete lattice. The various interactions 
between the groups are represented by ferromagnetic or antiferromagnetic couplings between nearest-neighbor spins 
(see Fig.|21l. Evidently, this is a very crude description of surfactant solutions and is not expected to yield quantitative 
predictions. Another difficulty is attaining thermodynamic equilibrium in simulations of these self-assembling systems, 
which contain slowly relaxing, large aggregates. Such models, however, have been shown to correctly reproduce various 
qualitative features of amphiphilic systems, e.g., aggregate formation, aggregate shape, and the overall structure of 
phase diagrams. The main advantage of this statistical-mechanical approach is that, by tuning a small number of 
parameters, one can get from the MC simulations insight into molecular mechanisms that determine the overall 
system behavior. Here we briefly present a model similar to that of Ref. |3l|. It will serve us when we discuss gemini 
surfactants in Sec. IIII Bl 
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FIG. 2; Schematic representation of a surfactant molecule (gray) solubilized in water (white) in a lattice spin model ;31]. Each 
water molecule and surfactant group are represented by a spin variable on a lattice site. Water molecules have spin -1-1, head 
groups +2, and tail groups —1. In between the head and the tail there is a "neutral" group of spin 0. The various particles 
interact via nearest-neighbor "ferromagnetic" couplings favoring spins of the same sign, except for the head-head interaction 
which is "antiferromagnetic" disfavoring neighboring heads. The chain connectivity and the overall number of chains are 
preserved during the MC simulation. 

In the lattice model each water molecule is assigned a spin tr = -f 1, a hydrocarbon group in the tail has spin cr = — 1, 
and the head group has spin a = +2 (see Fig. EJ. In between the hydrophilic head group and the hydrophobic tail 
there is a 'neutral' group of spin a = 0. All the couplings are ferromagnetic (favoring neighboring spins of the same 
sign), except for the head-head coupling which is antiferromagnetic, mimicking head-head repulsion due to screened 
electrostatics. The energy of the system can be written as 

i; = -J^(7,aj(l-2<5,^,2'5.,,2), (14) 

where J > is the coupling strength, (ij) denotes summation only over nearest-neighbor pairs of the lattice, and 
Sij = 1 when i = j and zero when z j is the Kronecker delta function. Two neighboring water molecules attract 
each other with energy — J because then UiGj = 1. The same applies to two tail groups, while a water molecule and 
a tail group repel with energy +J. The extra factor of (1 — 25ai,25aj,2) in Eq. l(T^ is unity for all cases except when 
i and j are two heads with a = -1-2, yielding a repulsion of -t-4J between two head groups. Finally, a head group 
and a water molecule attract with energy — 2J, and a head group and a tail group repel with energy -\-2J . These 
couplings apply whether the two neighboring lattice sites belong to the same molecule or not. In addition, the groups 
belonging to the same surfactant molecule are kept linked throughout the simulation. Thus, the essential features 
of hydrophobicity, hydrophilicity, molecular connectivity, and (screened) electrostatic repulsion are all accounted for 
using the single parameter J . Other parameters are the length of the tail group and total number of surfactant 
molecules in the system. 

The MC simulation starts from a certain configuration of surfactant molecules in water. At each iteration the 
various groups of the surfactants are moved while maintaining the connectivity of the molecules, their total number, 
and the total number of water molecules (e.g., using a "slithering snake" scheme HJ)- The energetic cost of the move 
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is calculated using Eq. H14() . and the MC step is accepted or rejected according to a Metropolis criterion, ensuring 
convergence toward equilibrium. 

This simple scheme can reproduce much of the richness of surfactant self-assembly, including the formation of 
monolayers, micelles and bilayers, the dependence of the cmc on tail length, transitions between various aggregate 
shapes, etc. On the other hand, such models can merely indicate general trends and not detailed information. For 
example, the correspondence between the MC spin variable representing sub-groups of a surfactant molecule and the 
actual chemical groups is not well defined and remains ambiguous to some extent. 

Another class of numerical studies that have been used to explore surfactant self-assembly are Molecular Dynamics 
(MD) simulations. These models range in detail from coarse-grained "bea d-sp ring" representations of the molecules 
(e.g., [Hill HIM mill) to atomistic descriptions (e.g., [S HI Ellil EJ ) • The advantage of the MD approach, 
as compared to phenomenological theories and spin models, is that the description of the system on the molecular 
scale is less artificial. The disadvantages are the limited spatial and temporal extent of the simulations, entailing 
equilibration problems, and sometimes also a large number of required parameters. A typical all-atom MD simulation 
of an aqueous surfactant system may contain about a hundred surfactant molecules along with few thousands water 
molecules, and the dynamics can be run for a few nanoseconds (e.g., Ref. f43|). A coarse- grai ned simulation allows 
a significant increase of these numbers at the expense of molecular detail (see, e.g., Ref. |33)- Here we outline a 
coarse-grained approach to surfactant micellization, as presented in Ref. |32l |. which was later extended to gemini 
surfactants 's^. 

The MD model of Ref. contains only two types of "particles": watcr-likc and oil- like, where a surfactant 
molecule is composed of a few water-like particles (the head group) and a chain of oil- like particles (the tail). The 
particles interact via a truncated Lennard- Jones potential. 



where r is the inter-particle distance, e the energy parameter of the Lennard-Jones potential, d its length parameter, 
and Tc a cutoff. This potential has a minimum at rmin = Hence, for rc < J^min the potential is purely repulsive, 

which is what is chosen for the oil-water interaction. For rc > rmin the potential contains a short-ranged repulsion 
followed by an attractive region, which is a suitable choice for the water-water and oil-oil interactions. In addition, 
the particles constituting a single surfactant molecule are connected by harmonic potentials of equilibrium length d 
and strong spring constant (much larger than e/rf^), ensuring chain connectivity. 

The MD simulation starts from a random distribution of surfactants in water. It then evolves in time according 
to the classical equations of motion governing the motion of individual particles. The simulations typically contain 
few ten thousands particles and are run for about 10^-10^ time steps |lllllE6, 37]. Thanks to the coarse-grained 
description, this can amount to about 1 /is in real system time |37|. Such a scheme was shown to successfully reproduce 
the structure of monolayers and micelles of various shapes, and to provide some understanding of the dynamics of 
surfactant self-assembly [H . On the other hand, as in the MC case, the coarse-grained representation prevents a 
well defined correspondence between the simulated system and the actual molecules in the experiments. 



C. Phase behavior 



Concentrated surfactant solutions and ternary water-oil-surfactant systems exhibit a rich variety of disordered and 
liquid-crystaUine phases 00,01. Some examples are the lamellar (La) phase, sponge (L3) phase, hexagonal {Hi) 
phase, and cubic (Vi) phase. All of these phases are based on various packing of surfactant layers — bilayers (in the 
binary-mixture case) or monolayers (in the ternary case). The lamellar phase is made of stacks of parallel layers, the 
sponge phase contains a disordered arrangement of multiconnected layers, the hexagonal phase consists of hexagonal 
arrays of parallel cylinders, and in the cubic phase the surfactant layers are spheres arranged in a cubic lattice. 

Unlike micellization one deals here with macroscopic bulk phases and their corresponding phase transitions. The 
powerful tools of thermodynamics and statistical mechanics are hence applicable. Consequently, the theory of surfac- 
tant phase behavior has reached a more advanced level, in particular, in the case of phases with long-ranged order. 
We shall not review these theories here, as most of them have not been used in current models of gemini surfactants, 
but merely mention the various approaches. 

Two phenomenological approaches to the phase behavior of surfactant binary and ternary mixtures have been used. 
The first is based on the Ginzburg-Landau formalism, which has been widely used in statistical physics | 4^ . It starts 
with a lattice description of the mixture and derives from it a coarse-grained, continuous expression for the energy 
(Hamiltonian) , which can be studied by statistical- mechanical techniques. The second approach is based on the elastic 
and thermodynamic properties of the membranes that make the various phases. For a review see, e.g., Ref. H^]. In 
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addition to these phenomenological theories, a variety of lattice spin models employing Monte Carlo simulations, as 
discussed in Sec. Ill B\ were originally designed and apphed to study surfactant phase behavior pa 1271 12M 1291 ISOl Isif . 



III. MODELS OF GEMINI SURFACTANTS 

Having provided the necessary background, we now turn to models of gemini surfactants. As will be demonstrated 
below, these models are essentially extensions of surfactant self-assembly theories, which have been reviewed in SecUll 
The binding of the surfactant molecules into pairs via spacer chains introduces new constraints affecting the molecular 
arrangement in monolayers, micelles and mesophases, as well as the thermodynamics of self-assembly. 

A. Surface behavior 

As in Sec.^ we begin by looking at a saturated monolayer, this time made of gemini surfactants, lying at the water- 
air interface. The gemini nature of the molecules, i.e., the introduction of the spacer, adds considerable complexity to 
the problem, mainly since it introduces anisotropy and inhomogeneity into the monolayer. A schematic view of the 
surface covered with dimers is shown in Fig. 13 The dimers may be oriented in various directions, and the distances 
between two linked monomers (in a dimer) and between two unlinked ones will differ in general. 



• • • 




• 




• • 





FIG. 3: Schematic top view of a saturated monolayer at the water-air interface. Left: regular surfactant head groups separated 
by a mean optimal distance ro. Right: Same view of gemini surfactants where the head groups are linked into dimers by 
spacers. The mean distance between head groups in a dimer is Vd, which in general differs from the distance between unlinked 
head groups. 



Nevertheless, we are going to disregard these complications and focus on the simplest question: how does the 
introduction of a spacer consisting of s groups affect the inplane distance r^{s) between two monomers belonging 
to the same dimer? Although seemingly oversimplified, the answer to this question will give us key insight into the 
surface behavior of gemini surfactants. 

We proceed by reviewing and slightly extending the work presented in Refs. |46L li^ . We have seen in Sec. Ill Al 
that the interaction between surfactants in a saturated monolayer can be roughly approximated by "effective springs" 
whose equilibrium length tq is determined by the optimum packing at saturation, and whose spring constant fco is 
given in Eq. Q . It is thus natural to consider the spacer chain as another "effective spring" , of equilibrium length 

and spring constant kg, connecting the two surfactant heads in a dimer. The combination of the two types of 
monomer-monomer interaction — the one present between unlinked monomers and the one due to the spacer — is 
then reduced to adding together two springs in parallel. From this we obtain 



fcpro -f- ks{s)rs{s) 



(16) 



The origin of the spacer "spring" is entropy, and its parameters are determined by the statistical distribution of 
spacer configurations. The "equilibrium length" of the spring is the mean end-to-end distance of the spacer chain, 
and the "spring constant" is inversely proportional to the variance of the end-to-end distance, 



kBT 
(r2) - (r)2 



(17) 
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where the averages are taken over all spacer chain configurations. Thus, the harmonic-spring approximation for the 
spacer is equivalent to representing the actual statistical distribution of spacer configurations by its first two moments. 

Before considering specific models for the spacer chain, let us examine what qualitative results are expected from 
this description. When the spacer is very short and rigid, such that fcg ^ fco, the equilibrium length of the dimer is 
determined by the spacer, ~ rg. On the other hand, when the spacer is very long and flexible, such that kg ^ ko, 
Td will be determined by the regular monomer-monomer interaction, rd — tq. Hence, upon increasing the number s 
of groups in the spacer, we expect rd(s) to first increase and then saturate toward tq, the optimal distance between 
the monomeric surfactants. Whether the behavior for intermediate spacer lengths is monotonous or not depends on 
specific details of the spacer chain. If the spacer stiffness fcs(s) drops sufficiently fast with s, the "interaction spring" 
will start dominating before rs{s) exceeds tq, and ^^(s) will then grow monotonously with s. By contrast, if fcs(s) 
decreases slowly with s, the "spacer spring" will dominate even for quite long spacers and — rg will become larger 
than ro- For even longer s, it will have to decrease back toward rg, leading to non-monotonous behavior in this case. 

The simplest model for the spacer is that of a Gaussian, constraint-free chain. This case is somewhat artificial and 
is discussed here merely as a model for very flexible and long chains, in contrast with more realistic models discussed 
below for more rigid chains. A Gaussian chain consisting of s segments is analogous to a random walk of s steps. The 
mean-squared displacement of such a walk, averaged over all s-step walks, should scale linearly with s. The mean 
end-to-end distance of a Gaussian spacer is therefore '--^ bs^^^, where b is the segment length. More specifically, the 
statistical distribution of the end-to-end distance in a Gaussian chain is 

p(r)dr = (^-rr) ^ e"^"-' /'■^"'''Unr^dr. (18) 



From the mean and variance of this distribution we get, according to Eq. H17|) 

isV^"^ . _ knT 



Thus, in order to calculate rg and fcg we just need to know the segment length b. For a polymethylene spacer b is 2.53 A. 
The remaining information required to compute rd from Eq. (|17|l are the properties of the monomeric surfactant in a 
saturated monolayer, namely fcp and rg. In a roughly hexagonal arrangement of molecules one has fco ~ 0.3 ksT/A^ 
(see Sec. Ill At . A saturated monolayer of DTAB surfactants, for example, is known to have ao ~ 55 A^, i.e., rg ~ SA. 
By using these values and substituting Eq. 1(1^)1 in Eq. rd as function of s is calculated and depicted as the 

dashed line in Fig.^ The inter-monomer distance increases moderately with s and even exceeds tq, yet the maximum 
and descent back to rp are shallow and occur at very large s, lying outside the experimentally relevant range of spacer 
lengths, 1 < s < 20 shown on the figure. 

More realistically, the spacer chain can be described by the rotational-isomeric model, where each segment in the 
chain can have only three possible orientations with respect to the two that precede it in the sequence (the three 
conformations called trans and gauche^) In addition, we require that the hydrophobic spacer be restricted to 

reside in the non-aqueous side of the water-air interface. These constraints stiffen the chain and bias its statistics 
toward larger end-to-end distances. We therefore expect a larger overshoot and a sharper maximum of rd(s), as is 
confirmed by the solid curve in Fig. 01 The points of this curve were obtained from simulations of rotational-isomeric 
polymethylene chains whose ends were fixed to a surface (th e air-water interface) and whose segments were forbidden 
from crossing that surface into the water side (see Ref. |4g for more details). From the simulations one obtains the 
end-to-end distance distribution and then extracts the "spring" parameters rg and kg according to Eq. l(T7|l |47l |. 

What has been calculated is the inter-monomer distance in a dimer and not the average area per dimer in the 
monolayer. However, because the latter must increase together with the former (cf. Fig.|3Jl, this very simple "spring" 
model reproduces the experimental observation of a non-monotonous behavior of a{s) for the m-s-m gemini surfactants 
and gemini surfactants derived from arginine Q. While the shape of the experimental curve is reproduced only 
Qualitatively, the position of the maximum at s ~ 12 is the same as the one found experimentally for 12-S-12 surfactants 
Furthermore, the spring model elucidates the source of the non-monotonous behavior — a competition between 
the regular monomer-monomer interactions on one side, and the natural length and rigidity of the spacer, on the 
other. According to this picture we should expect a more moderate and monotonous increase in a for more flexible 
spacer chain, as has been demonstrated by the Gaussian-chain example above. This may explain the behavior of a(s) 
observed for the m-EOa;-m gemini surfactants, having more flexible poly(ethylene oxide) spacers 4, 5]. (Compare, 
for example, our Fig. ^ with Fig. 1 of Chapter 4 in this volume.) 

These qualitative features, as well as the maximum at s ~ 10-12, were found to remain unchanged upon various 
refinements of the model, e.g., the inclusion of non-bonded interactions within the spacer chain, or a more detailed 
treatment of the monolayer structure |46j . A hydrophobic effect, i.e., repulsion of spacer monomers from the water 
phase, was invoked in several works as an explanation for a "lift-off" of the spacer from the water surface and hence 
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FIG. 4: The distance between the two head groups of a gcmini surfactant in a saturated monolayer as a function of the 
number of groups s in its spacer. The distance is rescaled by ro, the distance between unlinked heads in a saturated monolayer 
of the corresponding monomeric surfactant. Curves for two spacer models are shown: Gaussian chain with no constraints 
(dashed line) and rotational-isomeric chain restricted to the non-aqueous side of the interface (symbols and solid line). The 
latter has a maximum at s = 12, in accord with experiments 



the maximum in a{s). Such an effect, according to the spring model, actually suppresses the maximum as it brings the 
spacer ends closer together and thus reduces the overshoot of r^- We note that this effect might be related, though, 
to the maximum observed in the cmc of m-s-m surfactants at lower s values; see Sec. IIII Bl 

The description provided by the spring model is too simplistic to account for various details of gemini surfactant 
monolayers. In particular, two critical comments can be made. First, the experimentally observed decrease of a{s) for 
s > 12 is much steeper than what the model describes As has been suggested in Ref. 0, this might be a result of 
increased pre-micellar aggregation in the bulk solution as the spacer becomes more hydrophobic. Second, the model 
regards the spacer as an isolated chain, whereas in reality the gemini surfactant has two additional tail chains nearby. 
In this respect the model treats the geminis as equivalent to "bolaform" surfactants. While undoubtedly the presence 
of the tails is important for quantitative predictions, it is not expected to alter the qualitative competition picture 
described above, and a similar non-monotonous behavior of a{s) was indeed observed in bolaform surfactants as well 



B. Micelles 



The micellization behavior of gemini surfactants is qualitatively different from that of regular ones. We have 
reviewed some of these differences in the Introduction, and they can now be further elucidated in the light of what 
we have discussed in the previous sections. 

The cmc of gemini surfactants is typically one to two orders of magnitude lower than that of the corresponding 
monomeric surfactants 6] . The lower cmc can be directly attributed to the increase in the number of hydrocarbon 
groups in the molecule, i.e., decrease in the hydrophobic contribution g\^c, due to the second tail, and also due to the 
hydrophobic spacer chain in the case of m-s-m surfactants. Based only on the contribution of a second tail to g^c 
and the fact that the molecular volume is roughly doubled going from the monomeric surfactant to a gemini one, 
one would have predicted a larger decrease in cmc than what is actually observed. The difference is probably due to 
unfavorable terms introduced by the spacer, which will be further discussed below. 

The cmc of m-s-m gemini surfactants, instead of monotonously decreasing with the number s of spacer hydrocarbon 
groups (i.e., with molecular hydrophobicity) , is a non-monotonous function with a maximum around s ~ 4-6 [^0J3|. 
A corresponding non-monotonous behavior is observed in the Krafft temperature 9] and micellization enthalpy |lO|. 
This behavior can be attributed to the straightness and rigidity of short spacers, which force their hydrocarbon 
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groups to be in unfavorable contact with water. At about s ~ 4-6, although the spacer chain is still rigid, a gauche 
conformation should become accessible, allowing some of the groups to penetrate in the micellar hydrophobic core. 
When the spacer is hydrophilic this effect should be absent, as indeed is the case with m-EOx-m surfactants, exhibiting 
a weak monotonous increase of the cmc with the hydrophilic spacer length, x 

As a function of spacer length s, m-s-m surfactants exhibit an unusual progression of aggregate shapes from cylinders 
to spheres to bilayers. This is different from the more natural succession, occurring in monomeric surfactants, where 
the change in aggregate curvature is monotonous: spheres transforming into cylinders transforming into bilayers. 
Assuming that the molecular area at the aggregate surface is related to that in a saturated monolayer, this uncommon 
behavior can be qualitatively understood in view of the non-monotonous variation of o(s) as function of s, discussed 
in Sec. 1111 Al Considering that the radius and volume of the micellar core depend primarily on the tails and not on 
the spacers, an increase and then a decrease of a as a function of s should be accompanied by a decrease and then an 
increase in the packing parameter P of Eq. ©, hence the unusual morphological sequence. 

More specific predictions require a detailed theory and will be reviewed next. 



1. Phenomenological model 

An extension of the phenomenological theory of surfactant aggregation to gemini surfactants with hydrophobic 
spacers is presented in Ref. 50.J . It introduces the following additions and modifications to the model outlined in Sec. 
HTbI 

(i) The hydrophobic free energy (/he contains, apart from the double-tail contribution, also a spacer contribution. 
Only the spacer section which penetrates in the micellar hydrophobic core. Score, is considered. This spacer section is 
taken simply as the difference between the total spacer length and the mean head-head distance. Score — s —■ a}-^'^ jh. 
Since the three chains (two tails and spacer) are in partial contact already prior to aggregation, the hydrophobic 
energy gain per hydrocarbon group is taken to be smaller than in the case of a single-tail surfactant. 

(ii) The interfacial term, (?int(a)j is modified to account for the part of the core- water interfacial area that is 
now occupied by the spacer. The chain length that participates in this "shielding" is proportional to o}!"^ . This 
contribution is thus 

'Jffint - (72-7i)a'/'w^, (20) 

where 72 is the spacer-water interfacial tension and w is the spacer width. If the spacer is a polymethylene chain 
then 72 — 7i and this correction vanishes. 

(iii) When the spacer is short, it forces the two tails to be closer together than they would be if they belonged to 
two separate molecules. This packing constraint reduces the entropy of the tail chains. For a single-tail surfactant, 
the area close to the core-water surface sampled by tail groups is Otaii ~ vjR (cf. Fig.^l, with a prefactor of order 
unity that varies with aggregate shape |23l ISOj l . The proximity to a second tail due to the spacer reduces this available 
area per tail to Ogp (sfo)^. Thus, the contribution to the free energy can be estimated as 

(iv) Finally, the most difficult modification to handle is the electrostatic one. A short spacer forces the distance 
between two connected head groups to be shorter than that between two unconnected ones, resulting in a nonuniform 
charge distribution of pairs over the micellar surface (cf. Fig.|2Jl. This problem is bypassed in Ref. !50l| by introducing 
an empirical correction factor to gos, which becomes equal to unity when the spacer is longer than the mean inter-head 
distance. 

This extended phenomenological model is applied in Ref. 50] to gemini surfactants with short hydrophobic spacers, 
using parameters known from regular single-tail and double-tail surfactants. The model yields cmc values for various 
tail lengths in good agreement with the measured ones. More important, it correctly accounts for the observed micelle 
shapes of m-s-m surfactants with small s, i.e., the crossover from cylinders to spheres as s is increased. Following the 
changes in the various free energy contributions, one can identify the crossover mechanism as a competition between 
the elastic and packing contributions from the tails (favoring cylinders), and the electrostatic contribution (favoring 
spheres). Note, however, the various assumptions and approximations involved in these calculations. Although the 
electrostatic contribution to the free energy is found to be crucial for the self-assembly behavior, it is treated somewhat 
dubiously, as already admitted in Ref. [50|. 
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2. Computer simulations 

The additional complexity introduced by the spacers makes analytical calculations very difficult. One is inclined, 
therefore, to resort to computer simulations in order to gain detailed information on the self-assembly of gemini 
surfactants. 

The statistical-mechanical approach based on Monte Carlo (MC) simulations, as outlined in Sec. Ill Bl was extended 
to treat gemini surfactants |5lL l52l | . The "spin" assignment to various groups and the corresponding energy function 
are the same as for regular, monomeric surfactants [see Eq. (|14ll ]. The main modification is the connection of head 
groups in pairs via spacers (Fig. |SJl. Both hydrophobic spacers (spins a = —1) and hydrophihc ones (cr = +1) were 
simulated. In addition, the role of spacer stiffness was checked by assigning an energy penalty for "kinks" in the 
spacer configuration. 
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FIG. 5: Schematic representation of a gemini surfactant molecule (gray) solubilized in water (white) in a lattice spin model 
[51I l5^. The spin scheme is similar to that of Fig.|21 except that the two hydrophilic head groups (spin -1-2) of each surfactant 
are linked by a spacer chain. The spacer is composed of spins —1 for hydrophobic spacers, as shown here, or -1-1 for hydrophilic 
spacers. In addition, a "kink" in the spacer chain such as the one shown here is assigned an energy penalty in order to mimic 
the role of spacer stiffness. 

This spin model reproduces a few important properties of gemini surfactants as observed in experiments, primarily 
the non-monotonous dependence of the cmc on spacer length for hydrophobic spacers, and the formation of branched 
and entangled worm- like micelles in the case of short hydrophobic spacers (see Fig.EJ. However, the MC simulations 
produce also some findings which are not in full accord with experiments. The cmc is found to increase with tail length, 
unlike the common experimental results. (An exception to this rule is presented in the experiment of Ref. |53j . yet in 
a system having equilibration problems.) The mechanism for such a cmc increase with surfactant hydrophobicity is 
unclear. It is hard to simply attribute it to spacer-head repulsion, since the increase is found to be insensitive to the 
"spin" associated with the head group. A similar issue appears in the cmc dependence for hydrophilic spacers, which 
is found to decrease with spacer length, in disagreement with the experimentally observed (and expected) increase Q. 
The maximum in the cmc as a function of s is obtained for long hydrophobic spacers of about s ~ 12 regardless of tail 
length, contrary to the experimental result of only s ~ 5 Surfactants with long (s = 16) hydrophobic spacers 

are found to form rod- like cylindrical micelles, whereas in experiments they form bilayers !l4j. Hence, the spin model 
investigated by MC simulations seems to capture part of the essential features of gemini surfactant self-assembly while 
missing others. 

The "bead-spring" Molecular Dynamics (MD) approach discussed in Sec. lIIBl was extended as well to treat gemini 
surfactants j36l l37l| . The only essential modification is the connection of head groups in pairs by spacer chains. Like 
the tail chains, the spacers are made of "oil-like" particles connected to one another by harmonic springs, where only 
hydrophobic spacers were studied. These MD simulations are able to reproduce the micellar shapes formed by the 
m-s-m gemini surfactants — branched worm-like micelles and ring micelles — compared to the spherical morphology 
formed by the corresponding monomeric surfactants (see Fig. [TJ. A similar coarse-grained MD approach, along with 
a self-consistent-field calculation, were applied to the more complex glucitol amine gemini surfactants, which have 
flexible sugar side-chains attached to the charged head groups 54]. The main finding is a transition from cylindrical 
micelles to bilayers upon increasing pH, in accord with experimental indications 5Sj . 
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FIG. 6: Worm-like micelles formed in a Monte Carlo simulation of the spin model. The gemini surfactant has two heads of 
one lattice site each, two neutral groups of one site each, two tails of 15 sites each, and a hydrophobic spacer of two sites. 
Different gray tones correspond to different surfactant groups. The water molecules are not shown for clarity. (Reprinted with 
permission from Ref. |5lll.') 

C. Phase behavior 

The spacer length has an unusual effect also on the phase behavior of systems containing m-s-m gemini surfactants. 
This fact has been mentioned already in the Introduction. For binary surfactant-water mixtures, the regions in the 
concentration-temperature phase diagrams, where single-phase hexagonal and lamellar phases are the stable state, 
become smaller as s increases, vanish for s = 10-12, and then are finite again for s = 16 [T^ . In ternary water-oil- 
surfactant phase diagrams of the same surfactants, the size of the emulsification (single-phase) region increases with 
s and then decreases, with a maximum at s ~ 10 |l7l |. 

These observations can be rationalized in the light of the packing considerations discussed in Sees. IIIlXI and ITlI Bl 
As s increases from small values, the optimal area per molecule a{s) increases and the packing parameter P of Eq. 
© decreases. Hence, surfactant packing into bilayers, which are the building blocks of the lamellar and hexagonal 
mesophases, becomes less and less favorable. As we have seen in the previous sections, this effect is maximal for 
s = 10-12, whereupon bilayers apparently can no longer be stabilized. When s increases further, a decreases, and 
bilayers can form again. This qualitative description agrees also with the experimental results for m-EiOx-m gemini 
surfactants. The observed monotonous increase in the phase-diagram region belonging to the isotropic micellar phase 
|56| is in accord with the moderate monotonous increase of a{s) observed for these surfactants at the water-air 
interface Q. 

In ternary oil-water-surfactant mixtures, as s increases and P decreases, the surfactant monolayers required for 
stabilizing a microemulsion can have higher curvature. Hence, smaller oil domains can form and the microemulsion 
region in the phase diagram extends toward higher surfactant concentration. This effect, too, should be maximal for 
s = 10-12, as observed in the experiment 0. 

A theoretical study of gemini surfactant phase behavior, using MC simulations of a lattice model and a theory 
of mixture thermodynamics, is presented in Ref. 57] . Employing a simulation technique which, after equilibration, 
samples the composition of small regions in the entire lattice, the model is insensitive to long-ranged structures and is 
rather focused on the thermodynamics of phase coexistence. The study was restricted also to short hydrophilic spacers. 
Thus, the results cannot be compared with the experiments mentioned above. The main finding is the suppression 
of the three-phase region (coexistence of water-rich, oil-rich, and surfactant-rich phases) upon introducing molecular 
rigidity. 
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FIG. 7: Micelles formed in coarse-grained Molecular Dynamics simulations of gemini surfactants. Left: branched worm-like 
micelle. The surfactant has two heads of three water-like particles each, two tails of six oil-like particles each, and a spacer 
of one oil-like particle. Different gray tones correspond to different surfactant groups. The water molecules are not shown for 
clarity. (Reprinted with permission from Ref. [3^ . Copyright 1994 American Association for the Advancement of Science.) 
Right: ring micelles. The surfactant has two heads of one water-like particle each, two tails of four oil-like particles, and a 
spacer of two oil-like particles. Different gray tones correspond to different surfactant groups. The water molecules are not 
shown for clarity. (Reprinted with permission from Ref. |37|. Copyright 2002 American Chemical Society.) 



IV. CONCLUSIONS AND OPEN QUESTIONS 



The unusual self-assembly behavior of gemini surfactants poses challenging puzzles to theoretical investigations. 
We have reviewed the currently available models that attempt to address these puzzles, concentrating on surface 
properties, micellization, and phase behavior of gemini surfactant solutions. The overall impression emerging from 
the current state of the art is that, despite several successes, the theoretical understanding of gemini surfactants is 
fragmentary and lags behind the wealth of available experimental data. 

As demonstrated in this review, current gemini surfactant models are based on previous theories of surfactant 
self-assembly, with the most essential modifications required due to the addition of the spacer chains. It seems that 
this route has been exhausted, and further progress will depend on detailed consideration of features distinguishing 
gemini surfactants from regular monomeric ones. 

One of the distinct factors that stands out as a crucial ingredient is the spacer effect on lateral organization of 
the surfactant molecules in water-air monolayers and at aggregate surfaces. Linking the head groups in pairs has 
at least three different aspects as can be seen in Fig. O (i) The distribution of head groups on the surface becomes 
inhomogeneous as linked head groups have a mutual distance different from that of unlinked ones. This should affect, 
for example, the surface charge distribution, (ii) The spacers give the surfactant molecule an inplane orientation, i.e., 
the combined head group made of the two monomeric heads and spacer is very anisotropic. Such breakdown of isotropy, 
as is known from other systems, may lead to drastic effects on the overall behavior and may result in the formation 
of inplane liquid-crystalline order. Nematic ordering due to elongated head groups was theoretically addressed in the 
case of bilayer membranes ^58j . A recent work combining dichroism spectroscopy and atomistic MD simulations has 
revealed orientational ordering of gemini surfactants in cylindrical micelles 42]. More theoretical work is required to 
elucidate this issue, (iii) The above two aspects apply as well to double-tail surfactants (e.g., phospholipids) with 
large, elongated head groups. What truly distinguishes gemini surfactants from double-tailed surfactants is the fact 
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that the spacer makes a "soft" hnk between the head groups. Containing at least several chemical bonds, it allows 
a degree of conformational flexibility to the entire molecule, e.g., with respect to the relative orientations of the two 
tails. This feature is probably what allows gemini surfactants to form such uncommon structures as branched micelles 
and ring micelles. In the MD simulations of Ref. '36], for example, the gemini surfactants residing in a branching 
junction of a worm-like micelle were found to have their two tails oriented in different directions. This property might 
also make gemini surfactants serve as cross- linkers of regular micelles Hence, it looks like our understanding 

of gemini surfactant self-assembly will be incomplete until we have a good account of the interplay between various 
lateral organizations of these molecules at surfaces. 

Another important direction where there is substantial experimental information but almost no theory is the 
dynamics and rheology of gemini surfactant solutions. This aspect is particularly relevant to applications, as these 
solutions exhibit unusual and useful rheological properties such as shear thickening at low volume fractions P]. 
Moreover, recently these properties have made micellar solutions of gemini surfactants a model system for studying 
nonequilibrium behaviors such as shear thickening and ultra-slow relaxation [2H . We note that the dynamic issues 
and the issue of molecular organization mentioned above may be closely related. It has been argued recently that 
the distinct rheological behavior of worm- like micellar solutions (e.g., shear thickening) stems from the formation and 
interlinking of ring micelles 60]. 

We hope that this review and the posed open questions will motivate further theoretical studies of this class of 
fascinating and very useful self-assembling molecules. 
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